Electronic structure of underdoped cuprates 



Os 
OS 

S3 

a 

oo 
(N 



c3 



i 

Sh 

O 
o 



> 

so 
Os 



O 

Os 
03 



i 

Sh 
O 
O 



X 
S3 



Andrey V. Chubukov 1 ' 2 and Dirk K. Morr 1 

1 Department of Physics, University of Wisconsin, Madison, WI 53706 
2 P.L. Kapitza Institute for Physical Problems, Moscow, Russia 
(February 7, 2008) 

We consider a two-dimensional Fermi liquid coupled to low-energy commensurate spin fluctuations. 
At small coupling, the hole Fermi surface is large and centered around Q — We show that 

as the coupling increases, the shape of the quasiparticle Fermi surface and the spin-fermion vertex 
undergo a substantial evolution. At strong couplings, g 3> wo, where ujq is the upper cutoff in 
the spin susceptibility, the hole Fermi surface consists of small pockets centered at (±7r/2, ±7r/2). 
Simultaneously, the full spin-fermion vertex is much smaller than the bare one, and scales nearly 
linearly with \ q — Q\, where q is the momentum of the susceptibility. At intermediate couplings, there 
exist both, a large hole Fermi surface centered at (n, it), and four hole pockets, but the quasiparticle 
residue is small everywhere except for the pieces of the pockets which face the origin of the Brillouin 
zone. The relevance of these results for recent photoemission experiments in YBCO and Bi2212 
systems is discussed. 



I. INTRODUCTION 

The physics of cuprate superconductors has been a 
very popular issue for nearly a decade following the dis- 
covery of high temperature superconductivity by Bed- 
nortz and Miiller in 1986 jl). Over the past few years, 
it became increasingly clear to the "high-T c commu- 
nity" that the mechanism of superconductivity is di- 
rectly related to the unusual normal state properties 
of the cuprates, particularly in the underdoped regime. 
The 63 Cu spin-lattice relaxation rate and spin-echo de- 
cay rate [||J|], uniform susceptibility Q, in-plane and c- 
axis resistivity |5[|^], all demonstrate a temperature and 
doping dependence which is different from the predic- 
tions of the conventional Landau Fermi liquid theory. 
A behavior different from a conventional Fermi-liquid 
theory has also been observed in angle-resolved photoe- 
mission experiments [Q-pJ, transport measurements [ [To| , 
and optical experiments on the underdoped cuprates. 
Characterizing and explaining this behavior is one of the 
major challenges presently facing the high-T c commu- 
nity. A number of theoretical approaches to the cuprates 
have been proposed in recent years [^2[-[l7|] . One of the 
approaches to the cuprates, which we advocate in this 
paper, is based on the assumption that the physics of 
high-T c superconductors is governed by the close prox- 
imity of superconducting and antiferromagnetic regions. 
Indeed, parent compounds of cuprate superconductors 
(LaiCuOn, YBa 2 Cu 3 OQ, SV2CUO2CZ2) are antiferro- 
magnetic insulators. Upon hole doping, the antiferro- 
magnetic order rapidly disappears and the system even- 
tually becomes a superconductor. There is plenty of evi- 
dence that medium-range magnetic correlations are still 
present even at optimal doping (the one which yields the 
highest T c ). The most direct evidence came from neu- 
tron measurements for optimally doped La\- x Sr x CuO^ 
which reported the observation of propagating spin- waves 
at scales larger than J p8[. Additional evidence for 



strong spin fluctuations even at optimal doping comes 
from NMR and magnetic Raman experiments jl9],^0) . 

The key element in the magnetic approach to cuprates 
is the assumption that the dominant interaction between 
fermions is the exchange of spin fluctuations |17y . In 
this respect, the spin- fluctuation approach resembles a 
conventional BCS theory with the only difference that 
the intermediate bosonic mode is a magnon rather than 
a phonon. There is, however, one principal difference 
between the magnetic and the phonon mechanisms - in 
the static limit, the exchange by spin fluctuations yields 
an effective interaction which is positive (i.e. repulsive) 
for all q: 



ife+9,- 
k,-k 



(1) 



where g is a coupling constant, and x(q) > is the static 
spin susceptibility. In this situation, the BCS-type equa- 
tion for the gap, 
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cannot be satisfied by an s— wave type solution = 
const. On the other hand, in the vicinity of an anti- 
ferromagnetic region, the spin susceptibility is peaked 
at or near the antiferromagnetic ordering momentum 
Q = (71", tt) such that Eq.(^) relates A& and Ak+Q- One 
can then look for a solution which satisfies A& = — Ak+Q 
in which case the overall minus sign in Eq.(^) disap- 
pears. For the 2D square lattice and near half filling 
(where superconductivity has been observed), this extra 
condition on the gap implies that it changes sign twice 
as one goes along the Fermi surface, and vanishes along 
the Brillouin zone diagonals. In the language of group 
theory, such a pairing state possesses d x 2_ y 2 symmetry 
(Afc oc (cos k x — cos k y )). At present, there is a large body 
of evidence that this pairing state is the pairing state in 
cuprates - the most direct evidence came from tunnel- 
ing experiments on YBCO |2M. Notice, however, that 
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there still exist some experimental data which at present 
cannot be explained within the d— wave scenario [ p2| . It 
has been argued in the literature that there is an ad- 
mixture of the s— wave component in the gap, which is 
very small for underdoped materials, but increases as one 
moves into the overdoped regime p^,^4|. This admix- 
ture may be due to the fact that the materials contain 
orthorhombic distortions which mix the s and d x 2_ y 2 rep- 
resentations pq ). 

The effective interaction in Eq.(Q) has been used to 
compute the superconducting transition temperature in 
the Eliashberg formalism 26 27^. For g ~ Q.QeV inferred 



(0.7C) 



from resistivity measurements, and xil) inferred from 
NMR data, these calculations yield T c ~ 10 2 K which is 
of the same order as the actual value of T c . The ap- 
plication of the same approach to underdoped cuprates, 
however, yields a T c which steadily increases as the sys- 
tem approaches half-filling, simply because the spin sus- 
ceptibility becomes more and more peaked at the anti- 
ferromagnetic momentum. Experimental data, however, 
show that T c decreases when the system becomes more 
and more underdoped and eventually vanishes even be- 
fore the system becomes magnetically ordered. 

Another piece of evidence concerning underdoped 
cuprates comes from recent photoemission experiments 
which measure the electron spectral function A(k, u>) — 
(l/TT)ImG(k,uj) and can therefore locate the Fermi sur- 
face in momentum space. These experiments show that, 
at optimal doping, the hole Fermi surface is large, cen- 
tered around (71", tt), and encloses an area of about half 
the Brillouin zone, in consistency with Luttinger's the- 
orem pg}| . For underdoped cuprates, however, the data 
by the Stanford (?],|2£j and Wisconsin || groups indicate 
that the Fermi-surface crossing is present only for the 
momenta close to the zone diagonal (k x w iky). No 
crossing has been observed near (0, tt) and symmetry- 
related points. Instead, the spectral function in these 
regions of momentum space has a broad maximum at 
100 — 200rneT^. This feature can be interpreted as the 
quasiparticle peak in a strongly interacting Fermi-liquid 
in which a substantial portion of the spectral weight is 
shifted from the coherent peak into the incoherent contin- 
uum. Interpreted in this way, the data indicate that the 
Fermi surface crossing in the underdoped cuprates exists 
only for the directions close to the zone diagonal. This is 
consistent with the idea, first put forward by Shraiman 
and Siggia |30|, that the hole Fermi surface in heav- 
ily underdoped cuprates consists of small hole pockets 
centered around (7r/2, 7t/2) and symmetry-related points 
(see Fig. 0) . Whether this idea is fully consistent with the 
observations is still a subject of debate. If pockets exist, 
then one should observe two Fermi surface crossings along 
the zone diagonal. Kendziora et al. argued that they did 
observe two crossings in their measurements on ZQK su- 
perconductor, and therefore are able to reconstruct the 
whole pocket Q. The data by Shen et al. |nj for the 
most heavily underdoped superconductors also indicate 
that the maximum in the spectral function first shifts to 
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FIG. 1. A small Fermi surface consisting of a hole pocket 
centered around (tt/2, tt/2). 



lower frequencies as one moves from (0,0) to (tt/2, tt/2), 
then disappears before the (tt/2, tt/2) point is reached, 
and then reappears and moves towards higher frequencies 
as one moves from (tt/2, tt/2) to (tt,tt). They, however, 
argued that it is difficult to determine whether the reap- 
pearance of the peak between (tt/2, tt/2) and (71", tt) is an 
observation of the second side of the pocket, or whether 
its reappearance possesses a structural origin [ ^2[ . 

Our interpretation of the photoemission data parallels 
the one made by Kendziora et al. - we believe that there 
is substantial evidence that the Fermi surface in the un- 
derdoped cuprates does consist of small hole pockets. A 
natural question one may then ask is whether the forma- 
tion of pockets is related to the reduction of T c . We will 
argue in the paper that this is, in fact, the case. Specif- 
ically, we will argue that in the doping range where the 
Fermi surface evolves from a large one to a small one, 
vertex corrections substantially reduce the pairing inter- 
action. How precisely this evolution occurs is the sub- 
ject of the present paper. Notice that the evolution of 
the Fermi surface cannot be continuous since the large 
and small Fermi surfaces are centered around different 
points in momentum space. In fact, our results show 
that the Fermi surface evolution should necessary include 
a Lifshitz-type phase transition in which the topology of 
the Fermi surface changes. 

We conclude the introduction with the discussion of a 
microscopic model for the cuprates. There are numerous 
reasons to believe that the low-energy properties of the 
cuprates are quantitatively captured by the effective 2D 
one band Hubbard model I33|-|3T 
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(3) 



Here a is a spin index, = c\ a Ci. a , and Uj is a hop- 
ping integral which we assume to act between nearest 
and next-nearest neighbors (t and t' , respectively). By 
itself, this one-band model is a simplification, since in 
"first principles" calculations one would start with the 
three band model for the CuO^ unit. However, it is gen- 
erally accepted that at energies less than about 5eV the 
hybridization between Cu and O orbitals is rather strong, 
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and one can effectively describe the system by a single de- 
gree of freedom per C11O2 unit |36|. The on-site Coulomb 
repulsion U which we use in Eq.(|3|) is close to the actual 
charge-transfer gap between the unhybridized Cu and O 
bands which is about 2eV . 

Furthermore, the Hubbard model contains only elec- 
trons but no spins. Spin fluctuations appear in this model 
as collective modes of fermions. To obtain these modes, 
one has to sum up the RPA series in the particle-hole 
channel [^-^9). As a result of the summation, the prod- 
uct of the two fcrmionic Green's functions is replaced by 
a spin susceptibility whose poles correspond to spin fluc- 
tuation modes. The susceptibility thus obtained takes 
the general form 



x«(q,w) 



(4) 



1 - Uxii(q,u) 

in which Xu(q., u>) is the irreducible particle- hole suscep- 
tibility, and i = x,y, or z. Similarly, one can obtain 
the spin-fermion vertex (i.e., the coupling between two 
fermions and one spin fluctuation), by starting with the 
original four-fermion Hubbard interaction term and com- 
bining one incoming and one outgoing fermion into the 
RPA series. One can then introduce an effective model in 
which fermions and spins are considered as independent 
degrees of freedom coupled by 
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where cr, are the Pauli matrices, and the spin propagator 
(i.e. the Fourier transform of the retarded spin-spin cor- 
relation function) is the RPA susceptibility x{l,<^)- The 
coupling constant g is equal to U in the RPA approxima- 
tion. In semi-phenomenological theories, this coupling, 
however, is considered as some effective input parameter 
p6| . The argument here is that g can be substantially 
renormalized from U due to diagrams not included in the 
RPA series. 

The spin-fermion model is a convenient point of depar- 
ture if one intends to study the effects of spin fluctuations 
on the electronic spectrum. We point out, however, that 
this model neglects a direct fermion-fermion interaction. 
We will see later in this paper that one actually needs 
this interaction as a necessary term to restore the Ward 
identities. 

The spin-fermion model can be further simplified if 
one assumes some phenomenological form for the full 
spin susceptibility rather than computing it in the RPA 
approximation. This procedure can partly be justified. 
The point is that only the imaginary part of the ir- 
reducible particle-hole susceptibility comes from an in- 
tegration near the Fermi surface where the fermionic 
Green's function is known on general grounds. The com- 
putation of the real part of x(9;^)j on the other hand, 
involves an integration over regions far from the Fermi 
surface. In these regions, the actual fermionic propaga- 
tors can differ substantially from their values for the non- 
interacting fermions. This in turn implies that the RPA 



approximation may be insufficient, and a phenomenolog- 
ical form of x with the parameters deduced from the 
experimental data might be more appropriate. 

We now briefly discuss the phenomenological form of 
the susceptibility. One can argue quite generally that in 
the vicinity of the antifcrromagnctic phase, the full spin 
susceptibility behaves as [Eof 
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where xq = c*£ 2 j £ is the correlation length (measured in 
units of a lattice constant), c sw is the spin- wave velocity, 
and ui s f = Cg W /2£ 7 where 7 is a spin damping rate. The 
parameters in the phenomenological susceptibility can be 
inferred from comparison to NMR and neutron scatter- 
ing data, as it was demonstrated by Millis, Monien and 
Pines [[l0| . The choice of a is indeed arbitrary as it can be 
absorbed into the coupling constant. Furthermore, as we 
said above, u> s f can be computed in an arbitrary Fermi 
liquid provided that the spin-fermion coupling is not too 
strong. In fact, a comparison of the calculated (within 
the above framework) and measured u> s f yields an infor- 
mation about the value of the quasiparticle residue at the 
Fermi surface in optimally doped cuprates |4l| . 

In the further discussion, we will consider both Hub- 
bard and effective spin-fermion models. In the next sec- 
tion, we review the large U spin-density-wave approach 
for a magnetically ordered state at low doping, and show 
that in the presence of long-range order, the full vertex, 
g e ff(k), in fact vanishes at the antiferromagnetic mo- 
mentum k — Q. Simultaneously, the hole Fermi surface 
consists of four small pockets. We then argue that both 
of these features should survive when the system looses 
long-range order, and disappear only at much larger dop- 
ing concentrations when the system effecti vely looses its 
short-range magnetic order. In Sec [II B we study the 
effective model which reproduces the evolution from a 
large to a small Fermi surface as th e stre ngth of the 
spin-fermion coupling increases. In Sec [II D we compute 



the corrections to the spin-fermion vertex and show that 
these corrections substantially reduce the vertex for the 
range of couplings when the Fermi surface has a small 
area. Finally, in Sec. IV we state our conclusions and 



point to unresolved issues. 



II. ORDERED STATE 

We begin our discussion with the magnetically-ordered 
state. Our emphasis here will be to obtain the form of 
the spin-fermion interaction vertex, and the shape and 
the area of the quasiparticle Fermi surface. 

A straightforward way to study the ordered state is 
to apply a spin-density- wave (SDW) formalism p8| , p2| . 
Suppose that the system possesses a commensurate an- 
tiferromagnetic order in the ground state. Then the 
z— component of the spin-density operator, S(q) — 
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(1/2) E k c k+q,a^a,f3Ck,/3, has a nonzero expectation value 



(tt, tt) 



In the SDW approach, one uses 



at q = Q 

the relation (£ k A+q^a) = -<E* C l+Q,| c *a) = ( S *) 
to decouple the quartic term in Eq.(||). The truncated 
Hamiltonian is then diagonalized by a Bogoliubov trans- 
formation to new quasiparticle operators 



Ck,a = U k a k ,a + Vkbk,cr, 

Ck+Q,a = sgn(a)(ukbk,a - v k a k ,a), 
where the fermionic Bogolyubov coefficients are 



u k 
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Vk = 
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Here e k ' = 



(e k - e k+Q )/2, e k = -2t(cosk x + cosk y ) - 
At' cos k x cos k y is the single particle dispersion, and E k = 

\J (4 _) ) 2 + A2 ' wncrc A = U(S Z ). In the cuprates, 
"first-principle" calculations yield t ~ 0.3eV and t' ~ 
— 0.06eF p3[ . Note that only contains t, and, 

hence, the Bogoliubov coefficients do not depend on t' . 
The self-consistency condition for A has the simple form 
l/U = Et-l/Er, where the summation runs over all 
occupied states in the magnetic Brillouin zone. 

The diagonalization using (0) yields two bands of elec- 
tronic states (the conduction and valence bands) with a 
gap 2 A ~ U in the strong coupling limit (U ^ t,t'): 
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where E^ = E k ± e 



(+) 1+) 



= (e* + e k +Q)/2, and the 



prime restricts the summation to the magnetic Brillouin 
zone. We can further expand under the square root and 



obtain E 



J(cos k x 



cos ky) 2 =F 4i' cos k x cos k y , 



where J = At 2 /U is the magnetic exchange integral. 

It is also instructive to present the expression for the 
fermionic Green's function. It now has two poles at f2 = 
E k v {E = E — ii, where fi is the chemical potential), and 
vanishes at f2 = e.k+Q- Namely, we have 
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Consider now how the Fermi surface evolves as < S z > 
and, hence, A increases. Suppose that we fix the doping 
concentration at some small but finite level and increase 
A from zero to the strong coupling value A ~ U. For 
A = 0, the Fermi surface for holes has the form shown 
in Fig ga - it is centered at (tt, tt) and encloses the area 
which is slightly larger than half of the Brillouin zone, 
in accordance with the Luttinger theorem. Switching on 
an infinitesimally small A doubles the unit cell in real 
space. As a result, the electronic spectrum acquires an 
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FIG. 2. Schematic evolution of the Fermi surface with the 
SDW gap A in the ordered phase of the t — t' Hubbard model. 
We set t' = — 0.2t. For A = (a) the Fermi surface is large 
and centered around (tt,tt). For A = + (b) additional pieces 
of the Fermi surface appear, which are the images of the orig- 
inal Fermi surface shifted by (±7r, ±7r). For intermediate A 
(c) one observes hole and electron pockets around (n/2,n/2) 
and (0, 7r), respectively. For large A (d) the Fermi surface 
just consists of hole pockets centered around (it/2, it/2). The 
dashed line indicates the boundary of the magnetic Brillouin 
zone. 



extra branch which is just the image of the original dis- 
persion but shifted by Q (see Fig. ||). The same also 
happens with the Fermi surface - it acquires extra pieces 
because a Fermi surface crossing at k necessarily implies 
one at k + Q (Fig. . The subsequent evolution of the 
Fermi surface proceeds as is shown in Fig. f§i,d. The de- 
tails of this evolution are model dependent and therefore 
are not that relevant. An essential point is that, at large 
A, valence and conduction bands are well separated, and 
u k' v k ~ 1 /2, i.e., Z « 1/2 for all momenta. Since a sum- 
mation of the electronic states with quasiparticle residue 
Z — 1/2 over the full Brillouin zone gives the same result 
as a summation of the states with residue Z = 1 over half 
of the Brillouin zone, nearly all states in the valence band 
are occupied, except for a fraction, which is proportional 
to the doping concentration measured from half filling. 
This in turn means that the hole Fermi surface should 
be small and centered around the minimum of E k . For 
negative t' which we only consider here, this minimum 
is located at (7r/2,7r/2) @,|| (see Fig. |d). The area 
of this small, pocket-like Fermi surface can be obtained 
either by modifying Luttinger's arguments to account for 
the existence of both, two poles and a vanishing numera- 
tor in the Green's function, or by a direct computation of 
the particle density N/V. Performing this computation, 
we obtain 
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FIG. 3. The electronic spectrum along the (0, 0) — (tv, n) 
direction (a) and the (0,7r) — (n, 0) direction (b) for inter- 
mediate A (see Figj^c). The shaded area corresponds to the 
filled electronic states. Note the presence of electron pock- 
ets around (0, 7r) and (n,0). The dashed line in (a) is the 
quasiparticle spectrum for A = 0. 



where Shoie and S e i ec t r are the closed areas of hole pock- 
ets and doubly occupied electronic states (the latter ap- 
pear at intermediate stages of the Fermi-surface evolu- 
tion, see Figs. and ||b). At strong coupling (A ~ U ^> 
t), the doubly occupied electron states disappear and we 
have Shoie = x. 

We now proceed with the computations of the spin- 
fermion vertices. In the SDW formalism, the low-energy 
spin fluctuations (magnons) are collective modes in the 
transverse spin channel. Consequently, they correspond 
to the poles of the transverse spin susceptibility. The 
latter is given in the SDW theory by the series of ladder- 
diagrams. Consider first the strong-coupling limit at half- 
filling. Then each ladder can only consist of conduction 
and valence fermions. Since the unit cell is doubled due 
to the presence of the antiferromagnetic long range order, 
we have two susceptibilities — one with zero transferred 
momentum, and one with the momentum transfer Q = 
(tt, 7t). The explicit forms of these susceptibilities are |3g] 
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_ 1 + 79 

X ± (q,q + Q,uj) = -S 



uj — Q q + iS uj + Q q — i8 



UJ — fl q + iS UJ + Q q — l5 



Further, a sequence of bubble diagrams can be viewed 
as an effective interaction between two fermions medi- 
ated by the exchange of a spin-wave. The spin-wave 
propagators are i(Te g (t)eJ(0)) w = (fl q — uj — iS)^ 1 and 
i(Te q (t)e q (0)}uj = (Q q + ui — iS)^ 1 , where e q (e q ) are the 
boson creation (annihilation) operators, subindex uj im- 
plies Fourier transform, and the momentum q runs over 
the whole Brillouin zone. A simple experimentation then 
shows that the form of the two susceptibilities are repro- 
duced if one chooses the following Hamiltonian for the 
interaction between the original fcrmionic operators and 
the magnons (38): 



(12) 



Here il q = 4JS^/l — j q is the magnon frequency. 



H, 



el— mag 



u 



k q 



+ f h C k+q, a C k+QA e -q- e ?) 



V-/3- 



where 



1 (1 



V2 \ l +J q 



1/4 



1 (I 



79 



V2 \l-79 



1/4 



(13) 



(14) 



Performing now the above Bogolyubov transforma- 
tion, we obtain the Hamiltonian for the interaction be- 
tween the magnons and the conduction and valence band 
fermions 

Hel~m.ag = X] X! \_ a \ak a iP ,k+q(i\ ® aa{k , q) 
k q 

+ b Lk b *P,k+q4 $bb(k,q) + 4ak b W,k+<l e l ®ab(k,q) 
+ h \ak a iP,k+ q e\ ®ba(k,q) + H.c. 5 a -p- (15) 

The vertex functions are given by 

®aa,bb{k,q) = U ± (u k u k+q - v k v k+q yi q + (u k v k+q 

- V k U k+q ^T] q , 

®abM(k,q) = U (u k v k+q +v k u k+l ^j-q q =F (u k u k+q 

(16) 



VkVk+q I Vq 



In the strong coupling limit which we actually study now, 
they reduce to 



<S>aaMk,q) = [±(4 ' + 4+9)% + (< 
$ab,ba(k,q) = U [l] q TVq] ■ 



(") 



(17) 



We see that there are two types of vertices: $ a 6,ba 
which describes the interaction between conduction and 
valence fermions, and <fr a a,bb which involves either valence 
or conduction fermions. The first vertex is virtually not 
renormalized by Bogolyubov coefficients. However, at 
large A, it involves high-energy conduction fermions and 
therefore can be omitted. The second vertex is the one 
relevant for the pairing mechanism since slightly away 
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from half-filling, both valence fermions can simultane- 
ously be on the Fermi surface. We see that this vertex 
is substantially reduced at large coupling - it is of the 
order of the hopping integral rather than of the order of 
U. Moreover, it is easy to verify that <I>fc,& in fact vanishes 
as the magnon momentum approaches Q. For q close to 
Q, we have, expanding Eq.(^), &bb oc t\q — Ql 1 ^ 2 . The 
vanishing of at q = Q is a consequence of the Adler 
principle: in the ordered state, magnons are Goldstone 
bosons and their interaction with other degrees of free- 
dom should vanish at the ordering momentum to preserve 
the Goldstone mode to all orders in perturbation theory. 
Note in this regard that the vertex which includes con- 
duction and valence fermions should not be included into 
the corrections to the spin propagator. Indeed, this ver- 
tex is already used in the RPA series which yields the 
spin susceptibility presented in Eq. (|l2|) , and to include 
it into the corrections will just be double counting of 
the same diagrams. From this perspective, the fact that 
^ab,ba does not vanish for q = Q is not a violation of the 
Adler principle. 

We now construct the effective vertex for the spin- 
fermion model. Since the factors r\ q and fj q are derived 
from the spin susceptibility, we have to absorb them into 
S q . Actually, only r] q is relevant as it diverges at q = Q. 
Eliminating this factor in $fcb, and substituting U by the 
bare coupling constant g, we find that the effective spin- 
fermion interaction behaves as 

9eff=g {u k U k+q - V k Vk+q) OZ t\q - Q\. (18) 

This result was first obtained by Schrieffer We see 
that as long as one can keep the conduction fermions 
far away from the Fermi surface (which is the case at 
strong couplings), the effective pairing interaction be- 
tween fermions is substantially reduced by a vertex renor- 
malization |Q. Consider now the full static interaction 
between fermions mediated by spin fluctuations. It is 
given by T(q) = (g ef f{q)) 2 xio)- Since X (q) oc (q - Q)~ 2 
and g e ff(q) oc \q — Q\, T(q) tends to a positive constant 
as q approaches Q. On the other hand, however, we 
would need a substantial enhancement of the interaction 
near Q for the appearance of d— wave superconductiv- 
ity. We see, therefore, that the presence of a long-range 
order actually excludes the occurrence of d— wave super- 
conductivity at large U. In the next section we will study 
in detail what happens when the sublattice magnetiza- 
tion decreases and the mean- field gap, A = U < S z >, 
becomes smaller than the fermionic bandwidth despite a 
large U. At this point, however, we merely emphasize 
the correlation between the relevance of vertex correc- 
tions and the shape of the Fermi surface, namely, that 
when the Fermi surface consists of small hole pockets, 
vertex corrections arc relevant and the full vertex is much 
smaller than the bare one. 



FIG. 4. The lowest order self-energy corrections for the va- 
lence fermions. The solid and dashed lines are the bare prop- 
agators for the conduction and valence fermions, respectively. 
The wavy line describes transverse spin fluctuations. 

A. Fluctuation effects 

Our next goal is to study how fluctuations modify this 
simple mean-field picture. We will show that the mean- 
field form of the quasiparticle dispersion remains quali- 
tatively (and even partly quantitatively) correct even at 
strong couplings though there is a substantial shift of 
the spectral weight from the coherent part of the spec- 
tral function to the incoherent part which stretches upto 
frequencies of a few t. At the same time, the relation 
A = U < S z > does not survive beyond mean-field the- 
ory, and the strong-coupling, SDW-like form of the quasi- 
particle dispersion with a gap between two bands holds 
even when the system looses long-range antiferromag- 
netic order. More specifically, we will show that there ex- 
ist two different self-energy corrections. One self-energy 
correction is strong but nonsingular. It gives rise to a 
shift of the spectral weight towards the incoherent part of 
the spectral function, but preserves coherent excitations 
on the scale of J. The second correction renormalizes 
the gap but does not change the form of the quasiparti- 
cle Green's function. This correction is however singular 
for vanishing < S z > and breaks the relation between A 
and < S z > such that the gap survives and remains of 
0(U) when < S z > vanishes. To perform a perturbative 
expansion around the mean-field SDW state, we need an 
expansion parameter. Obviously, an expansion in U will 
not work as we assume that U 3> t. There exists, how- 
ever, a formal way to make the mean-field theory exact - 
one should extend the original Hubbard model to a large 
number of orbitals at a given site, n c — 25, and perform 
an expansion in 1/5 The 5 = oo limit corre- 

sponds to the mean-field solution, and all corrections are 
in powers of 1/5. The mean-field result for the spectrum 
is the same as in Eq.(^) apart from the extra factor of 25 
(we also have to redefine i — t/2S as t now scales with 
5). 

Consider now the lowest-order self-energy correction to 
the propagator of the valence fermions. To first order in 
1/5, there are two self-energy diagrams, one involves only 
valence fermions, and the other involves both, valence 
and conduction fermions (Fig. ||). In the first diagram, 
the vertex is reduced from U to 0(t) and vanishes at 
q = Q. Numerically, however, this diagram yields large 
corrections since the large term, A ~ US, in the energy 
denominator is cancelled out because both fermions are 
in the valence band. The denominator is then of the order 
of O(JS). As a result, the self-energy correction from 
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this diagram scales as 0(P/JS) — 0(U) and obviously 
contains some dispersion at this scale [H. At the same 
time, the bare dispersion is at the scale of JS. We see 
that one clearly needs more than just smallncss of 1/5 to 
make perturbation theory work - the ratio U / JS should 
be small too. This last assumption is not justified for 
high-T c materials where (U/JS) is typically of the order 
of - 20 - 30. 

We recently considered this diagram in detail and per- 
formed self-consistent calculations of the electronic dis- 
persion assuming that 1/Scl but also U/JS » 1 [ |49| . 
One can show that in this limit one can neglect vertex 
corrections but has to include the whole scries of rain- 
bow diagrams for self-energy corrections [Q. This cor- 
responds to the insertion of the full Green's function into 
the internal fermionic line in Fig^. The self-consistency 
equation for G(fc,<x>) then takes the following form: 



MF 



where 



MF 



d 2 q^{k iq )G(k + q,LU + n q ) (19) 



A + 2 JS(cos k x + cos k y ) 2 — At' cos k x cos k y 



is the mean-field dispersion and 



= 32PS v 2 k + vl +q - 2v k v q v k+q 



+ \ 1 - V q{ V k+q - V k ) I J 1 - V f 



Note that Eq. ( pL 9| ) is an integral equation in both momen- 
tum and frequency space. 

The form of the quasiparticle dispersion in the ordered 
phase is interesting in its own because of recent pho- 
toemission data on Sr^CuOiCl^ ]5l]j52] |. These data 
show that (i) the valence band dispersion is isotropic 
near the top of the band which is at (tt/2,tt/2), and 
relatively well described by the mean-field formula with 
if = — 0.5J, but (ii) coherent excitations exist only near 
(ir/2, ir/2). As soon as one deviates from this point, the 
width of the quasiparticle peak rapidly increases, and at 
the scale of 2 J, the excitations become mostly incoherent. 
Monte-Carlo and finite cluster simulations in the insulat- 
ing phase have similarly demonstrated that there exist 
coherent excitations at scales smaller than 2J, but also 
incoherent excitations which stretch upto a few t [^3|,[54| . 

Our analytical results are consistent with the experi- 
mental data and with the numerical simulations. We first 
discuss the critical value of ui for which one first observes 
a finite imaginary part of G(k,ui), or equivalently, the 
onset of a finite density of states (DOS) pq] . On a mean- 
field level, this happens on the scale of oj + A = O(JS). 
The analysis of Eq.(|l9]) shows, however, that the onset 
of a DOS actually occurs at the much larger scale of 
to + A = 0(t^/S). At these energy scales, the frequency 
shift due to fl q on the r.h.s. of Eq.(|ig|) can be neglected, 




FIG. 5. Schematic form of the total DOS resulting from 
the self-consistent solution of Eq.[H| W is the bandwidth for 
free fermions. 



and we just have to solve an integral equation in momen- 
tum space. Doing this by standard means, we obtain 
that a nonzero DOS appears at 



uj cr + A = ±4.2iV25. 

In the vicinity of ui cr , the DOS behaves as N(oj) ~ 
y/\u> — uj cr \. We see, therefore, that the incoherent part 
of the Green's function stretches upto scales comparable 
to the original fermionic bandwidth. 
Another important issue is the form of the quasiparticle 
Green's function at I UJ cr — CO | < JS, where one can show 
that the contribution from the magnon dispersion in fact 
cannot be neglected. In addition, photoemission experi- 
ments have shown strong evidence for coherent fermionic 
excitations in this region. We therefore assumed the fol- 
lowing trial form of the fermionic Green's function in the 
vicinity of ui cr 



G(k,u) = 



U! — U! cr + A(k) — i^{u! — LU cr ) 2 



(20) 



Performing self-consistent computations in this region, 
we found p9| that the quasiparticle residue is indeed 
small, Z ~ 0{Jy/S/t) < 1. Analog ous result has been 
obtained by Kane, Lee and Read [p0[. At the same time, 
A(k) is zero right at (tt/2,tt/2) and scales quadratically 
with the deviation from this point. The scale for the 
coherent dispersion is JS - the same as for the mean- 
field solution. Even more, we recovered the isotropic dis- 
persion near (ir/2,ir /2), namely A(k) = 2JSk 2 , k — 
(ir/2, 7r/2) — k, and obtained the same bandwidth 2 J as 
in the experimental data. For the damping rate we found 
7 ~ 1/JS. It follows from our results that there exists 
just one typical energy scale JS which separates coher- 
ent excitations at smaller energies and incoherent exci- 
tations at higher energies (see Fig||). It is essential that 
this holds in the strong-coupling limit t 3> JyS when 
the quasiparticle residue Z<1, and the mean-field dis- 
persion is completely overshadowed by self-energy cor- 
rections. In this respect, the existence of coherent exci- 
tations at the scales of JS turns out to be an intrinsic 
property of the Hubbard model and not the result of 
a mean-field approximation. Moreover, since the vertex 
does not depend on t' , the strong-coupling results are in- 
dependent of the t'/J ratio, contrary to the mean-field 
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dispersion [|56[ . Notice also that the dominant corrections 
due to valence- valence vertex come from short-range spin 
fluctuations. The contribution from long-range fluctua- 
tions with q « Q is reduced due to a vanishing vertex 
at q = Q. In other words, the valence- valence diagrams 
are certainly non-singular. Not surprisingly then, these 
corrections can reduce the coherent part of the quasi- 
particle dispersion but cannot qualitatively change the 
mean-field description. For the sake of simplicity, we will 
just neglect these corrections for the rest of the paper. 

We will now discuss the self-energy correction which 
involves both conduction and valence fermions. Here the 
vertex is of the order of U , but the energy denominator 
scales as US. As a result, one obtains a momentum- 
independent correction to the gap, and also (doing ex- 
pansion in J/U and in u>/U, where w is the external fre- 
quency) a correction to the dispersion on the scale of 
O(J). The latter is nothing but the conventional 1/5 
correction to the mean-field dispersion which is obviously 
small at large S. In other words, the inclusion of the self- 
energy diagram with valence and conduction fermions 
leaves the structure of the quasiparticle Green's func- 
tion unchanged, but renormalizes the gap. It is essen- 
tial, however, that the gap renormalization can be made 
singular when < S z > vanishes. The easiest way to see 
this is to consider what happens at infinitesimally small 
temperatures in 2D. The Mermin- Wagner theorem tells 
us that immediately as the temperature becomes finite, 
the long-range order should disappear. In perturbation 
theory, the onset of this effect is the appearance of log- 
arithmically divergent correction which scales as T log L 
where L is the system size. We evaluated the self-energy 
correction at T = + and obtained 



a) 



N ^ 



1 - 



1 - v 2 



1 - V 2 

q 



(l + 2n b q ) 



(21) 



where n b q = (exp(Sl q /T) — 1) _1 is the Bose distribution 
function. At any finite T, E is logarithmically divergent: 
E cx T J d 2 q/q 2 ~ TlogL. If this were the only cor- 
rection, the mean-field description would be completely 
destroyed. However, there exist another logarithmically 
divergent diagram which cancels the contribution from E. 
The point is that the mean-field expression for the gap 
Amf = U < S z > contains the exact sublattice magne- 
tization. On the mean-field level, one has (at large U), 
< S z >« S. When T becomes nonzero, < S z > also 
acquires a logarithmically divergent correction which in- 
volves the valence-conduction vertex (see Fig. |qJd) . In 
analytical form, we have 



(l+2n b q ) (22) 




b) 




Combining these two terms, we find that the logarith- 
mically divergent corrections to A are cancelled out. As 



FIG. 6. The loop diagrams for < S z >. The diagram in (a) 
represents < S z > /S on the mean-field level while those in 
(b) correspond to the lowest order corrections to < S z > / S 
due to the exchange of transverse fluctuations. 



a result, A remains finite even when the long-range or- 
der disappears. Moreover, a comparison of Eq.(21) and 
Eq. ( p2f ) shows that the gap remains exactly equal to US 
in the large U limit. 

By itself, this result is not surprising since at large 
U , the valence and conduction bands are nothing but 
the lower and upper Hubbard bands, respectively. In 
the near atomic limit, the gap between the two bands 
should be equal to U (for S — 1/2) upto corrections of 
0((t/U) 2 ) independent of whether or not the system is 
ordered. We extended our calculations to analyze the 
form of the full quasiparticle Green's function and found 
that not only A but also Gf u u(k,u>) is free from logarith- 
mical singularities and therefore does not undergo sharp 
changes when one looses long-range order due to singular 
thermal fluctuations. In other words, we found that the 
SDW structure of the electronic states with two coherent 
bands separated by a large gap rs U survives when the 
system looses long-range order. We remind that the key 
features of the SDW state at large U are the small Fermi 
surface and a strong reduction of the pairing vertex. We 
see that both features may exist even without a long- 
range order. The interesting issue then arises how the 
system evolves with doping in the paramagnetic state and 
how it eventually restores a large Fermi surface and a mo- 
mentum independent vertex function. Notice that with- 
out strong magnetic fluctuations, the two bands which 
emerge from the Hubbard levels are likely to be mostly 
incoherent, and there also appears, at finite t, a coher- 
ent band at about U/2. This was found in the infinite- .D 
studies of the Hubbard-model J57J . Before we proceed to 
the discussion of the Fermi surface evolution, we consider 
another, complimentary approach to the ordered phase, 
introduced by Kampf and Schrieffer ]5§|] . They demon- 
strated that one can obtain a mean-field SDW solution 
in the spin-fermion model without introducing a conden- 
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FIG. 7. (a) The Gorkov equation for the normal Green's 
function for the SDW state. The single and double solid lines 
correspond to the bare and full fermionic propagators, (b) 
The self-energy diagram which includes two anomalous loops 
The straight dashed line represents the interaction term U. 
(c) Formal representation of the same diagram as in (b) with 
the exchange of a 8— functional longitudinal spin fluctuations. 



sate, but by assuming that the longitudinal spin suscep- 
tibility has a 8— functional peak at zero frequency and 
momentum transfer Q: Xiil^) = (1/4) 5(u)5(q — Q). 
To recover the mean-field SDW result within the spin- 
fermion model one has to compute the lowest-order self- 
energy diagram which includes only the exchange of lon- 
gitudinal fluctuations. This diagram is shown in Fig.J^c. 
We then obtain 



E(k, w ) = x 

4 uj — e 



k+Q 



(23) 



where, as before, e k = e k — /i. Substituting this diagram 
into the Dyson equation, one immediately recovers the 
mean-field SDW result for the spectrum with g being 
the equivalent of U. 

Furthermore, one can also reproduce the correct form 
of the vertex between fermions and transverse spin fluctu- 
ations. The bare vertex is equal to g. To restore the cor- 
rect mean-field form of the full vertex, one has to include 
the second-order vertex renormalization due to longitudi- 
nal spin fluctuations. The second order vertex correction 
is given by the diagram in Fig ^. The evaluation of this 
diagram is straightforward, and we obtain 



6g 



l 



4 (w - e k+Q ) (uj + fl - e k+Q+q ) 



(24) 



where q is the magnon momentum, and u> and f2 are the 
external fermionic and magnon frequencies, respectively. 





FIG. 8. The diagrammatic representation for the full ver- 
tex (filled triangle) between transverse spin fluctuations and 
fermions in the ordered state. The solid and dashed wavy 
lines correspond to the transverse and longitudinal spin fluc- 
tuations, respectively. The solid straight lines are the bare 
fermion Green's functions. Note that the momenta of the 
incoming and outgoing fermions differ by Q — (tv, n). 



Combining this result with g, we find the total vertex in 
the form 

(u> - e k+Q )(uj + Q- e k + q +Q) - (g/2) 2 
9tot — 9 t —rj Tri = \ \ zb ) 



(u> - e k ){uj + Q - e k+q+ Q) 



For the vertex which involves only valence fermions, we 
have at the mass surface, ui = —E%, uj + fl = —E^ +q . 
Substituting this into Eq. ( p5| ) and doing elementary ma- 
nipulations, we obtain the same expression as in Ec 
This result has also been obtained by Schrieffer 

One may wonder what happens with the higher-order 
self-energy and vertex corrections. In fact, they simply do 
not exist. The reason is that the use of the 8— functional 
form of the longitudinal susceptibility is just a way to re- 
express the mean-field decoupling without formally intro- 
ducing the condensate. Specifically, one can rewrite the 
Gorkov equations for the SDW state by introducing the 
condensate, but without formally introducing the anoma- 
lous Green's function, as is shown in Fig. [7]. The self- 
energy diagram obtained in this way is shown in Fig. 0b. 
We can now formally combine the two anomalous loops 
into the 8— functional longitudinal susceptibility in which 
case we obtain the diagram in Fig. 0c. However, this is 
just a formal way to reexpress the standard mean-field 
results. Still, any inclusion of the renormalization of the 
internal fermionic line in fact means that one inserts more 
pairs of anomalous loops into the self-energy, which one 
cannot do since this will render S reduccablc. 

The Kampf-Schrieffer approach has the advantage over 
the conventional SDW decoupling in that it can straight- 
forwardly be extended to the region where long-range 
order is not present, but where strong antiferromagnetic 
fluctuations exist, i.e. where the correlation length is 
large, and the spin susceptibility is strongly peaked, but 
not divergent, at zero frequency and at the antiferromag- 
netic momentum. These are the basic assumptions for a 
nearly antiferromagnetic Fermi-liquid which we are going 
to study in the following section. 







a) 




a) 



b) 



t 



•♦■ i i 4* i i 



■IT 

t : t : t 



b) 



o 




FIG. 9. faj The lowest-order self-energy diagram for the 
spin-fermion model without long-range order. The wavy line 
describes the exchange of spin fluctuations. All three fluc- 
tuation modes contribute equally to the self-energy, (b) The 
self-energy diagram in the FLEX approximation. The double 
line describes the full fermionic propagator. 



III. THE NEARLY ANTIFERROMAGNETIC 
FERMI-LIQUID 

A. Transverse and longitudinal modes 

At first glance, the calculations performed in the or- 
dered phase using the spin-fermion model can be directly 
extended to the case of a "nearly" 5-functional form of 
the spin susceptibility in the disordered phase. However, 
this extension requires some care since in the disordered 
phase, one no longer can distinguish between transverse 
and longitudinal fluctuations. Both modes behave in the 
same way, and both diverge as one approaches the anti- 
ferromagnetic transition. 

Now, if we apply the second-order diagram to evaluate 
the Green's function (see Fig.^,), we do obtain the SDW 
form of the electronic spectrum as Kampf and Schri- 
effer have demonstrated, provided that the coupling is 
larger than the upper cutoff, ojq, in the spin susceptibil- 
ity. However, because all three modes (two transverse 
modes and one longitudinal mode) equally contribute to 
the self-energy, the second-order self-energy at large cou- 
pling g 3> ujq is given by 

3 

H a {k,u)=g 2 22<r^a a 0a / Xaaiq^G^k + q^ + n) 

a=l 

«^ 5 2 G (fc + Q,^) (26) 

(we normalized J Xaa(q,ty to 1/4). This self-energy 
yields a "near" SDW solution with the relative correc- 
tions of the order (cu /g) 2 , but with the gap larger by 
a factor of -\/3 than the one in the ordered phase (i.e., 
A 2 w 3(g/2) 2 ). At the same time, the vertex correction 
term at g 3> does not acquire an extra overall factor: 
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FIG. 10. The counter diagram which cancels the contribu- 
tion from the lowest order self-energy term due to transverse 
spin fluctuations in the ordered state. The transformation 
from a to c is explained in the text. The solid and dashed 
lines represent the bare fermionic propagator and the inter- 
action U, respectively. The solid and dashed wavy line de- 
scribe transverse and longitudinal spin fluctuations, respec- 
tively. Observe that this diagram still contains the direct 
fermion-fermion interaction U. 



3 

Sg/g = g 2 22 a ^i a ^ a s0 / G Q {k+p,uj + Q) 

a=l •* 

x G Q (k + q +p, lu + Q + to) Xaa(p,u) dpdu) 
m -(g/2) 2 G (k + Q,oj) G Q (k + q + Q,uj + n) (27) 

As a result, the relative vertex correction on the mass 
shell (w « - A, G ~ 1/A) will scale as Sg « -g (g/2A) 2 
and will only account for a reduction of g by a factor 
1/3. This in turn implies that the full vertex will not 
vanish for q — Q, contrary to what one should expect 
if the mean-field SDW solution survives the loss of the 
sublattice magnetization. The reason for this discrep- 
ancy is that we have not yet included the counterterm 
which in the ordered phase cancelled out the correction 
due to transverse fluctuations thus eliminating two out 
of the three components of the spin susceptibility in the 
self-energy diagram. Indeed, if we leave just one compo- 
nent of x in Eq-(H^)i we find for the gap A w (g/2) in 
which case Sg rs —g and one recovers the correct form of 
the vertex. 

In order to identify the diagram which eliminates two 
components of the spin susceptibility in the paramag- 
netic phase, we go back to the ordered state, reformulate 
the mean-field SDW theory as a result of the exchange 
of spin fluctuations and find the diagram which cancels 
the contribution from the transverse susceptibility. Ob- 
viously, this counterdiagram should include a correction 
to < S z >. The corresponding diagram for the Hubbard 
model is shown in Fig. |io|a. 

Note that the internal lines in the anomalous loop are 
the full quasiparticle Green's functions given by Eq. (|o|) . 
Performing the same manipulations that lead to the ex- 
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change of longitudinal fluctuations (i.e., connecting the 
two anomalous bubbles by U (Fig.|lO|b) and then substi- 
tuting the single bubble by the full RPA series) , we arrive 
at the diagram shown in Fig. [Io| c. This diagram includes 
spin-fcrmion vertices for the exchange of both, transverse 
and longitudinal fluctuations, but also the bare Hubbard 
U. The presence of the bare U is crucial - it implies 
that this diagram cannot be obtained in the pure spin- 
fermion model since it neglects a direct fermion-fermion 
interaction. One can easily check that for a (5-functional 
form of the longitudinal susceptibility this diagram and 
the second-order diagram with transverse spin exchange 
cancel each other in the large U limit leaving the diagram 
with longitudinal spin exchange as the only relevant self- 
energy diagram. 

Suppose now that the long-range order is lost but spin 
fluctuations are strong. Then our reasoning is the follow- 
ing. Assume first that only one component of the sus- 
ceptibility should be counted in the self-energy. Then, 
as we just discussed, we obtain an SDW solution for 
g = U » u>o- We then substitute the SDW form of 
G(k,ui) into the remaining part of the self-energy (the 
same as in FigJ^a, but with only two components of the 
susceptibility left) and into the fermionic bubble of the 
diagram in Fig. |lO| c. Computing the total contribution 
from these two diagrams we found that this contribution 
is small compared to the one from Fig. f^c by a factor 
(wo/ff) 2 - This in turn justifies the use of only one compo- 
nent of the susceptibility in the second-order self-energy. 
Note that at small q the situation is different: the coun- 
terdiagram in Fig. [LOp is of higher order in g and can 
be neglected. In this limit, we recover a conventional re- 
sult of the lowest-order paramagnon theory that all three 
modes of spin fluctuations contribute to the self-energy. 

We now briefly discuss higher order corrections. In the 
paramagnetic phase, they indeed exist and, moreover, 
are not small at large g. However, extending the same 
line of reasoning as above to higher-order diagrams, one 
can easily verify that the higher-order self-energy and 
vertex corrections nearly cancel each other in the large g 
limit. For example, the third-order self-energy and vertex 
correction terms are shown in Fig. O. Assuming that just 
one fluctuation mode contributes to the renormalization 
of the internal fermionic line, we find after simple algebra 
that the two corrections nearly (to order (uio/g) 2 ) cancel 
each other. One, therefore can restrict with only the 
second-order diagram and hence fully recovers the mean- 
field SDW solution. 

Before we conclude this section we want to comment on 
what happens if we restrict our calculations to the pure 
spin-fcrmion model. In this case each inclusion of the self- 
energy correction introduces a factor of 3 due to the ab- 
sence of a counterterm. At the same time, the inclusion of 
the vertex correction only yields a relative factor of one. 
This implies that the self-energy corrections are more rel- 
evant than the corrections to the vertex. This point can 
be made rigorous by extending the SU (2) spins to SU (N) 
and taking the limit of large N |59). For general N, one 
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FIG. 11. Third order diagrams for (a) the self-energy and 
(b) the vertex correction. At large g these two diagram nearly 
cancel each other. 



substitutes the Pauli matrices a in (§) by the (A 2 - 1) 
traceless generators of SU(N). In the large N limit, the 
vertex corrections are small by a factor of 1/(N 2 — 1) 
and their relative contribution disappears at N — > oo. 
At the same time, higher-order self-energy corrections do 
not contain extra powers of 1/N compared to the second- 
order diagram, and therefore should all be included. The 
full self-energy correction is then given by a series of rain- 
bow diagrams, which obviously can be reexpressed as a 
second-order diagram with the full Green's function for 
the intermediate fermion (see Fig||b). This approxima- 
tion is known as the fluctuation exchange (FLEX) ap- 
proximation (the fully self-consistent FLEX approxima- 
tion also uses the full fermionic Green's functions in the 
RPA series for the spin susceptibility). The solution of 
the FLEX equations always yields a large, Luttinger-type 
Fermi surface with progressively decreasing quasiparticle 
residue as the spin-fermion coupling increases. We argue 
that at least in the large coupling limit, this procedure is 
incomplete because for each inclusion of the rainbow dia- 
gram one should also include a counterterm which effec- 
tively eliminates N 2 — 2 components of the susceptibility 
out of the self-energy correction. Performing calculations 
along this line, one indeed recovers the SDW results. 



B. The model 

We now discuss the model we are going to study. In a 
"first principles" calculation, we would have to consider 
the Hubbard model. At large U, this model contains 
two peaks in the density of states (the upper and lower 
Hubbard band) for any doping concentration. Far away 
from half filling, however, the total spectral weight of the 
upper band is small, and the excitations in this band 
are likely to be incoherent. As the system approaches 
half-filling, the spectral weight is shifted from the lower 
to the upper band and one gradually recovers the SDW 
form of the spectrum even before the system becomes 
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magnetically ordered. How this evolution occurs is one 
of the key issues in understanding the normal state of 
the cuprates. We have not yet solved this problem in the 
Hubbard model, instead, we considered the evolution of 
the spectral function in a toy spin-fermion model. We 
assume that the density of holes is fixed at some small 
but finite level, and vary the coupling constant g. In 
doing this, we in fact model the system's behavior as 
it approaches half-filling simply because by all accounts, 
the spin-fermion coupling should increase as the system 
becomes more "magnetic" . 

Further, for any g, we compute the self-energy cor- 
rections restricting with just one component of the spin 
susceptibility. The argument here is that, at strong cou- 
pling, a counterterm which we discussed in the previous 
section, cancels out the corrections due to the other two 
components of the susceptibility. At small couplings, we 
indeed will miss the overall factor of 3 in the self-energy, 
but this does not seem relevant as we do not expect any 
qualitative changes in the fcrmionic spectrum. On the 
other hand, the study of the toy model should give us the 
answers to two key questions: (i) how does the Fermi sur- 
face evolve from a large one, centered at (jr,ir), at small 
couplings to a small one, centered at (7r/2,7r/2) at large 
g, and (ii) how does the spin-fermion vertex evolve from 
a momentum-independent one at small g to one which 
vanishes at Q (upto (ui /g) 2 ) terms) at large couplings. 

Finally, we assume that the susceptibility has the form 
ofEq.(§) with some given oj s f , £ and c sw which do not de- 
pend on the strength of the spin-fermion coupling. This is 
indeed an approximation. When g is large and precursors 
of the SDW state are already formed, the particle-hole 
bubbles which constitute the RPA series for the suscep- 
tibility should indeed be computed with the full Green's 
functions and the full vertices. However, as we will see 
below, at large g, we in fact probe the scales larger than 
luq, in which case the structure of x(q, u/) at u> < ujq 
is irrelevant - we only use the fact that the susceptibil- 
ity obeys the sum rule. At weak coupling, g < u>q, the 
form of the susceptibility is relevant, but in this limit the 
real part of the susceptibility is just an input function, 
independent on g, while the imaginary part of \ can be 
computed using the bare fermionic Green's functions and 
vertices which, in fact, we will do later. 

The form of \u (Qi w ) as m Eq. (|^) is indeed an expan- 
sion near u> — and q = Q. Meanwhile, xil^ 1 - ) should 
indeed satisfy the sum rule J duj d 2 c[xu(q, w) = 1/4. To 
impose the sum rule, one should either add extra terms 
into Eq.(^) with higher powers of (q — Q) 2 , or impose a 
cutoff in the momentum integration. For computational 
purposes, it is easier to impose a cutoff luq. We have 
checked that the results for the Fermi-surface evolution 
practically do not depend on whether we impose a cutoff 
only in the momentum integration or also in the inte- 
gration over frequency. In the latter case, however, the 
computations are much easier to perform and a number 
of results can be obtained analytically. Below we present 
the results for the model with a cutoff luq in both mo- 



mentum and frequency. 

After presenting the model, we now proceed to our cal- 
culations. We will compute the full quasiparticle Green's 
function and the full vertex restricting with the second- 
order diagrams only. Here we apply the same reasoning 
as before, namely that in both the weak and the strong 
coupling limit, the relative corrections to the second- 
order expressions are small and scale as (g/ujg) 2 and 
(wo/ff) 2 , respectively. This implies that the second-order 
theory will yield a correct limiting behavior at weak and 
strong couplings. At intermediate couplings, g ~ luq, 
higher-order terms are indeed not small. However, even 
without them, we found in our numerical studies that for 
intermediate g, the solutions for the Green's function and 
the vertex strongly depend on the values of cusf,^, etc., 
i.e. the behavior is highly nonuniversal. It is therefore 
unlikely that the inclusion of higher-order diagrams will 
give rise to new behavior which is not already present in 
the second-order theory. 

We deem it essential to point out that, though we 
restrict the corrections to the second-order diagrams 
only, the computations of G f u u(k,u>) = Go(fc,w)/(l — 
£(£;, Lv)Go(k, uj)) (Go is the fermionic Green's function 
at g = 0) and the full vertex contain a nonperturba- 
tive self-consistency procedure which is necessary to de- 
scribe a transformation from a large to a small Fermi 
surface. Namely, we consider the chemical potential fi 
as an input parameter in the calculation of the self- 
energy, and then determine it from the condition that 
the total density of particles is equal to a given number, 
J du;d 2 kGf u u(k,u}) = (1 — x)/2. The nonperturbative 
nature of this procedure appears since the self-energy is 
not small at large g, and one finds a finite region in [k, uj) 
space where T,(k, uj)Go(k, w) > 1 i.e. the perturbative ex- 
pansion of G fuii(k,u)) in powers of £ does not converge. 
This region, on the other hand, contributes to the den- 
sity of particles, and this makes the evaluation of /a a 
nonperturbative procedure. 

We are now ready to present our results. We begin 
with the self-energy corrections which give rise to the 
Fermi surface evolution. Some of the results presented 
in the next subsection have already been reported earlier 



C. Fermi-surface evolution 

As before, we assume that without a spin-fermion in- 
teraction, the system behaves as a Fermi-liquid with 
a large Fermi surface which crosses the magnetic Bril- 
louin zone boundary (see Fig.|l2]a and 13 1). Near the 
Fermi surface, the bare fermionic Green's function is 
Go(k,uj) = Z i 1 (uj — (efc — //)), where £k — — 2t(cosk x + 
cos k y ) — At' cos k x cos k y . The chemical potential at g = 
is clearly of the order of t': /.i = At' cos 2 k J where fc£ is the 
x— component of the momentum for the point where the 
Fermi surface crosses the magnetic Brillouin zone bound- 
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FIG. 12. The evolution of the hot spots with the coupling 
strength. For g < g£) (a) there exist only one hot spot; for 
g = gcr' it splits into three (b) . As g continues to increase, one 
progressively looses two of the hot spots (c) , and for g ^> g^, 
(d) only a single "new" hot spot is present 



ary. These points are generally known as "hot spots" J6l| . 
In the Matsubara formalism, which is more convenient for 
computations, the self-energy at T = is given by 



S(fc,w m ) = J d 2 qdtt rn x{q,tln 

x G (k + 



(28) 



As we discussed above, we will use the phenomenologi- 
cal form of the susceptibility in Eq.(||) with a cutoff ojq 
in both momentum and frequency. The location of the 
Fermi surface is obtained from 

GjZ u (k,u m =0) = -e fc -£(fc,w m = 0) = (29) 

It is instructive to consider first the case when the spin 
damping is absent, because then we can obtain a full an- 
alytical solution for E. Since the susceptibility is peaked 
at Q, we can expand the energy of the internal fermion as 
f-k+q = e k+o +Vk+Q{q — Q). Substituting this expansion 
into Eq.([28|) and integrating over momentum and fre- 
quency, we obtain a rather complex function of the ratio 
Q>k = Vk+Q/c sw which we present in the Appendix. The 
formula for u) m ) simplifies considerably if a% = 1. 
Then we have 
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FIG. 13. Fermi surface evolution with increasing g. We 
used (in units of t), A = 0.1, ujq = 0.3, t' = -0.45, x = 0.1. 



For these parameters we obtain g 



(i) 



,( 2 ) 



1.64. 
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The figures are for g = 0, 2uo > S > ffc 

(2 ) ( 2 ) ( 2 ) 

g > Jcr , g > gcr , and g 3> g C r , respectively. In (c), the 
nested pieces of the Fermi surface are shown in bold. In Fig. 
(e), we also presented the values of the quasiparticle residue 
along the Fermi surface. Notice that the definition of g used 
in this paper differs by a factor of 2 from the one in Ref . M] . 



E(fc,w ri 



4 



1 



~k+Q 



iu> m — ek+Q 



(31) 



if 



k+C 



< Wn. Here we introduced Aq 



We have checked numerically that the qualitative pic- 
ture for the evolution of the Fermi surface does not de- 
pend substantially on the value of a& . We will therefore 
discuss the Fermi surface evolution by restricting with 
the form of the self-ener gy g iven in Eqs.(|30|) and (J3l|). 

The self-energy in Eq. ( p0|) is precisely what we need to 
recover the SDW solution which, we recall, yields a small 
Fermi surface centered at (ir/2, n/2). However, to obtain 
this small Fermi surface, Eq.(|3(]) should be satisfied (at 
zero frequency) for all points in fc— space, and, in par- 
ticular, for the points along the magnetic Brillouin zone 
boundary. At g — 0, the hot spots are at Ik — £fe+Q = 0. 
Clearly, at small g, the location of the Fermi surface 
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near the hot spots is obtained with the self-energy from 
Eq. (|3ll) . A simple experimentation shows that the loca- 
tion of the Fermi surface crossing does not change with 
small enough g, i.e., it still occurs at Ek — £fe+Q = 0. 
However, as g increases, there appear two new hot spots 
at efe ^ 0. To see this, consider a point at the mag- 
netic zone boundary right near a hot spot. Expanding 
Eq.© in e k +Q and substituting the result into Eq.fl29|), 
we obtain 



Gf ull {k,Q) = -e k 



9 



8A 



o(v / ^" 



A 2 - A ' 



(32) 



For sufficiently small g, the only solution for the Fermi 
surface is Ik = 0. However, when g reaches the value 
g£) = 2w [2A /( v /wfTAf + A )] 1/2 , the velocity along 
the magnetic zone boundary vanishes. For larger g, one 
still has a Fermi-surface crossing at ik = 0, however, two 
new hot spots appear (see Fig.[l3|b) with 

( (gilY 



1/2 



y/uj% + A 2 , - A 



1 - 



(33) 



As g increases, |e/c+Q| for the new hot spots also in- 
creases, and at g — 2ujq, these new hot spots sat- 
isfy |efe+g| = t^Oj i- e -j they become the solutions of 
G fuii( k > 0) = with the SDW-like form of the self-energy 
in Eq.([j0|) rather than the form in Eq.©. These solu- 
tions have a very simple form: e k = ig/2. For even 
larger g, the SDW form of the solution for the Fermi sur- 
face extends beyond the hot spots and is located between 
points D and D' in Fig.fl3"|c. As g continues to increase, 
the D and D' points from neighboring hot spots approach 
each other. Finally, at 



9 = 9. 



< 2 > 2cj 



8t(|/*| -u ) 



1 1/2 



2t+^ 2 -4|i'|(H-u;o) 



these points merge and the system undergoes a topo- 
logical, Lifshitz-type phase transition (Fig.|l3|d) in which 
the singly-connected hole Fermi surface splits into hole 
pockets centered at (±7r/2, ±7r/2) and a large hole Fermi 
surface centered around (ir, n). As g increases further, 
the large Fermi surface shrinks (Figjl3p) and eventually 
disappears at which point the Fermi surface just consists 
of four hole pockets (Fig.|l3|f). 

We found numerically that the topological phase tran- 

(2) 

sition at g cr is accompanied by drastic changes in the 
functional behavior of the system. To demonstrate this 
point wc plot in Fig. |l4|a,b the chemical potential /i and 
the area of the electron states in the Brillouin zone as a 
function of g. One clearly observes that both quantities 

are practically constant upto g cr , but increase consider- 

(2) 

ably for g > g cr (we attribute the small variations be- 
low gi? in Fig. |l4"| b to numerical errors). In other words, 
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FIG. 14. (a) The chemical potential as a function of the 
coupling, g. The parameters are the same as in Fig. |l^. For 

(J. = fio 



free fermions, 



T.16. The arrow indicates the value 



of gcr when hole pockets are formed, (b) The area of the 
occupied states vs g. The dashed line is the area for free 
fermions. 



(2) 

Luttinger's theorem is satisfied below g cr , but is violated 
above this critical coupling. We attribute the violation 
of Luttinger's theorem to the nonconvergence of the per- 

(2) 

turbative expansion in g above gcr . We remind in this 
regard that, in essence, Luttinger's proof is perturbative: 
he showed that I = J Gf u udH/du = order by order 
in perturbation theory. Implicit in this proof is the as- 
sumption that the perturbative series is convergent. We 
computed I numerically and found that it is equal to 

( 2) 

zero (within the accuracy of our calculations) below g cr 
but rapidly increases as soon as g exceeds the critical 
value. We reserve a detailed discussion of this issue for a 
separate publication. 

Further, we mentioned in the introduction that it is 
still a subject of controversy whether the sides of the 
hole pockets which are facing (tt, it) have been observed 
in experiments or not. The relevant physical quantity 
here is the quasiparticle residue at the Fermi surface. 
For small g, when one has a large Fermi surface, Zk is 
almost k— independent and close to one. For very large 
g, when one effectively recovers the SDW form of the 
electronic excitations, Zk is again k— independent and is 
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equal to 1/2. For intermediate g , however, Zj. is strongly 
k— dependent and anisotropic. The full expression for 
Zk is presented in the Appendix. The results for the 
case when hole pockets and a large Fermi surface coexist 
are presented (for a& = 1) in Fig.^3^. We see that the 
quasiparticle residue along the large Fermi surface is very 
small, which makes the experimental observation practi- 
cally impossible. Furthermore, we see that the quasipar- 
ticle residue of the side of the hole pocket which is facing 
(•7T, 7r) is roughly five times smaller than the one in the 
momentum region facing the T point. This may very well 
explain the experimental difficulties in observing the part 
of the hole pocket facing (it, ir). 

Finally, we want to discuss how the above results 
change when we take a damping of spin fluctuations 
7 = c 2 w /(2u> s f^ 2 ) into consideration. We found that the 
general scenario of the Fermi surface evolution does not 
change, but that its onset occurs at substantially larger 
values of the coupling constant than in the absence of 
damping. Specifically, we found (still, assuming for sim- 
plicity that a,k = 1) that the value of gir where a hot 

spot splits into three is gil\-f) — gir (7 = 0)^(7/Ao) 
where *(x) = 1 - x/(3n) + 0(x 2 ) for x < 1, and 
^{x) = (16 \ogx/ (nx)) 1 ' 2 for x 1. In the nearly an- 
tiferromagnetic Fermi liquid, 7/A0 = c su ,/(2w s /£) is a 

large parameter. In this case, g^ cx [7/ log^/Ao)] 1 / 2 
which is substantially larger than g$ oc (Ao) 1 / 2 which 
we obtained in the absence of damping. 



D. Vertex corrections 

The bare interaction vertex in the spin-fermion model 
is a momentum-independent constant g. At small cou- 
plings, the full vertex is indeed close to the bare one. At 
strong couplings, on the other hand, we have precursors 
of the SDW state, and, as we discussed, the full vertex 
should be much smaller than the bare one. We now study 
how the full vertex evolves with the coupling strength. 
Contrary to naive expectations, we find that the evolu- 
tion of the vertex is not smooth. For dcfinitcncss, we 
will study below the vertex which describes the interac- 
tion between fermions at the hot spots and transverse 
spin fluctuations with antiferromagnetic momentum Q 
and zero frequency. This vertex is relevant to the pairing 
problem since the incoming and outgoing fermions can 
simultaneously be placed on the Fermi surface. 

Consider first the limit of weak coupling. The diagram 
for the vertex correction is presented in Fig. ||. Since the 
spin susceptibility is peaked at Q, the internal fermions 
are also located near the hot spots, and all three inter- 
nal lines in the diagram carry low-energy excitations. A 
simple dimensional analysis shows that the vertex cor- 
rection is logarithmically singular in the limit when the 
external fermionic frequencies are zero and the correla- 
tion length is large: Ag oc g * log[maa:(a;i ! 2, c st0 ^ -1 )]. 
In order to obtain an analytical expression for Ag, we 
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FIG. 15. The graphical representation of the angle (f>o be- 
tween the normals to the Fermi surface at hot spots (dashed 
lines). For clarification we omitted the parts of the Fermi 
surface in the second and fourth quadrant. 



expand the fermionic energies near the hot spots as 
£fc = v h (k- k h ) cos(0) and e k +Q = Vh(k — kh)cos(4> + (f>o) 
and perform the integration over frequency and momen- 
tum in the spin susceptibility. Keeping both, the loga- 
rithmically divergent and the £- independent terms in Sg, 
we obtain for zero external fermionic frequency 



Ag 
9 



g 2 z 2 XQU S f 




log[sin(0/2)] 

COS (j) + COS 4>o 

0{8 2 ) (34) 



where S = c 2 sw /{2^ v h £) = £ui s f/v h < 1. Z is the 
fermionic quasiparticle residue in the absence of the spin- 
fermion interaction and (f>o is the angle between the nor- 
mals to the Fermi surface at k and k+Q (see Fig. |l5|). For 
the t — t' model for the quasi partic le dispersion we fin d 
0o = tt/2 + 2tan- 1 [{\ - 2\t'\ ju~JW)/{{\ + 2\t'\ ^JW)]. 
Observe that Ag does not depend on the upper cutoff in 
the frequency and momentum integration (i.e., on loq). 

We see that Ag diverges logarithmically when £ di- 
verges and S oc £ tends to zero. Alternatively, one can 
compute Ag at = but at finite tt>i,2, and also ob- 
tain a logarithmical divergence. This last result was also 
obtained by Altshuler et al. fl5|| . 

The integral in the right hand side of Eq. (|34|) can eas- 
ily be evaluated numerically (and also analytically for 
particular values of 4>o)- It turns out that this integral 
is negative for all </>o, i.e., the relative vertex correction 
is positive. This implies that for small couplings, ver- 
tex corrections in fact increase the coupling strength and 
therefore enhance the d— wave pairing interaction. This is 
the opposite behavior of what we would expect at strong 
couplings. The strength of the vertex correction is an- 
other issue which recently was the subject of some con- 
troversy |l^,^,^]. The key issue here is whether one 
should consider u s f as an independent input parameter, 
or assume that the damping is mainly due to the interac- 
tion with fermions. In the latter case, which we believe is 
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closer to reality, u;~^ by itself is proportional to g 2 Z 2 jv\ 
in which case the coupling strength, fermionic velocity, 
and the bare quasiparticle residue are cancelled out in 
the r.h.s. of Eq.(|34|). In explicit form we found for iv s f 



TT 

4 



uj sf = — | sin 00 1 -n 



9 2 z 2 xq ' 



(35) 



Substituting ui s f into Eq.([34|) we obtain, neglecting terms 
ofO(<5 2 ) 



sm 0o | 
4tt 2 



Re 



log[sin(0/2)] sin(0/2) 
log ■ 



COS + COS 0o 



6 2 



(36) 



This expression only depends on 0o. Numerically, the 
r.h.s. of Eq.(^) is small for all realistic values of S. 
Thus, for S = 0.27, t' = -0.45*, and \i = -1.16 (which 
corresponds to 10% doping), we have A o/q » +0.04. A 
similar result has been obtained in Ref. |Q . 

Notice that though the vertex correction is small nu- 
merically, Eq. ( |36] ) shows that the perturbative expansion 
for Ag breaks down since Ag exhibits a step-like behav- 
ior when g becomes finite. This behavior of Ag, how- 
ever, is just an artifact of our approximation in which 
we neglected the u 2 term compared to Lu/uj s f in the spin 
susceptibility. A more detailed analysis shows that a con- 
tinuous behavior of Ag is restored for g < 5. 

Consider now what happens when g increases and the 
Fermi surface starts to evolve. As we discussed in the pre- 
vious subsection, the evolution begins at g = gir with 
the flattening of the Fermi surface at the hot spots. For 
g > gir there appear two satellites of the original hot 
spot (see Fig. |l2|b). The central hot spot is a solution of 

ek = fJ, below and above gir , and the vertex corrections 
for this hot spot are virtually insensitive to the onset 
of the Fermi surface evolution (we recall that we restrict 
ourselves with the second-order correction which involves 
only bare fermionic propagators). On the other hand, for 
the two new hot spots, — /i is finite (see Eq.(33)), and 
thus the vertex corrections will be different. As the Fermi 
surface continues to evolve, the old hot spot and one of 
the two new hot spots eventually disappear leaving a sin- 
gle new hot spot as the only remaining one on the Fermi 
surface. The vertex correction at this hot spot depends 
on the value of the chemical potential and changes sign 
when \fi\, which grows with g for g > g„ , becomes com- 
parable to the fermionic bandwidth. At even larger g, 
\/j.\ grows approximately as g/2. In this limit, the rela- 
tive vertex correction becomes, to leading order in too/g, 
^9/9 ~ ~(.9/m) 2 J d 2 qdflxaa(q,^) ~ -1, i-e., the total 
vertex g to t = .9 + Ag is nearly zero. Performing calcu- 
lations beyond the leading order in LUo/g and for spin 
momentum q ^ Q, we found that in this limit, 



gtot(q) oc — 

g 



A 2 f^) +(q-Q) 2 



1/2 



(37) 



where A = O(l) is a numerical factor. We see that for 
very large couplings, the vertex is nearly linear in (q — Q) 
as it should be in the SDW phase. 

The strong coupling form of the vertex is very similar 
to the one suggested by Schrieffer E(| , the only difference 
is that in his expression, (ujo/g) 2 in Eq.(^7|) is replaced 
by £~ 2 . We argue that one needs more than just a large 
correlation length to recover the precursors of the SDW 
state in the electronic dispersion and the vertex. Namely, 
in our approach, the SDW form of the vertex appears due 
to a separation of the original fermionic dispersion into 
two subbands separated by an energy scale g which is 
large compared to the cutoff loq in the spin susceptibility. 
When the correlation length becomes larger, ojo clearly 
goes down (and hence, g/oJo increases), but luq remains 
finite when £ becomes infinite. When the system becomes 
ordered, there appears a new momentum scale, go o c < 
S z >, such that when \q — Q\ < qa, gtot(q) is strictly 
linear in \q — Q\ in accordance with the Adler principle, 
while for larger momentum, one obtains a crossover to 
Eq. (|37|) . For q » (tj /g), the difference between the 
two limits becomes negligible and one fully recovers the 
strong-coupling SDW result for the vertex. 



IV. CONCLUSIONS 

We first summarize our results. The goal of this paper 
was to study the evolution of the electronic and magnetic 
properties of cuprates as one moves from optimal doping 
into the underdoped regime. Specifically, we were mostly 
interested in a possible relation between the reduction of 
T c with decreasing doping and the loss of the pieces of the 
Fermi surface near (0, tt) and related points, as evidenced 
in photoemission data. To address this issue, we consid- 
ered an effective spin-fcrmion model in which itinerant 
fermions interact with low-energy spin fluctuations whose 
dynamics are described by a semi-phenomenological spin 
susceptibility which we assume to be peaked at the an- 
tiferromagnetic momentum. In the absence of the spin- 
fermion interaction, the fermions are assumed to form a 
Fermi liquid with a large Fermi surface which encloses an 
area roughly equal to half of the Brillouin zone. We ar- 
gued that one can model the system's evolution towards 
half-filling by increasing the strength of the spin-fermion 
interaction. The weak coupling limit models the situa- 
tion near optimal doping, while the strong coupling limit 
corresponds to strongly underdoped cuprates. We found 
that, as the coupling increases, the Fermi surface first 
evolves in a continuous way, and its area remains un- 
changed. At some critical coupling, however, the system 
undergoes a topological, Lifshitz-typc phase transition in 
which the singly connected hole Fermi surface centered 
at (tt, tt) splits into a hole pocket centered at (tt/2, tt/2), 
and a large hole Fermi surface centered at (tt, tt). As 
g increases further, the area of the large Fermi surface 
gradually decreases, and it eventually disappears. As a 
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result, the Fermi surface at large enough g consists of 
four small pockets. Simultaneously with the changes in 
the Fermi surface topology, the spin-fermion vertex also 
changes from being nearly constant at small g to being 
nearly a linear function of q — Q at strong couplings. This 
last form of the vertex is the one expected for an ordered 
antifcrromagnet. When the vertex is linear in q — Q, one 
does not obtain an attractive interaction in the d— wave 
pairing channel since the enhancement of the spin sus- 
ceptibility near Q is fully compensated by the reduction 
of the vertex. 

Our results, therefore, show that there is an inter- 
play between the Fermi surface evolution and the loss 
of d— wave superconductivity in strongly underdoped 
cuprates. Both phenomena are related to the fact that 
the strong coupling SDW forms of the electronic disper- 
sion and the spin-fermion vertex do not change dras- 
tically when the system looses long-range magnetic or- 
der, contrary to what one could expect from the mean- 
field theory. Instead, the electronic structure gradually 
changes towards a conventional Fermi liquid as magnetic 
correlations become less and less relevant. We believe 
that the Fermi liquid picture of electronic states, with 
one single peak in the density of states, is restored some- 
where around optimal doping. 

As we already mentioned in the introduction, the 
transformation from a large to a small Fermi surface in 
underdoped cuprates is consistent with the experimental 
data, though only one group so far reported the obser- 
vation of both sides of the hole pocket. We emphasize, 
however, that in the intermediate coupling regime, which 
most probably corresponds to the experimental situation 
for T c = 30K and T c = 60K Bi22Yl superconductors 
studied by photoemission , the quasiparticle residue 
for the part of the pocket which faces (it, tt) is rather 
small (few times smaller than on the other side of the 
pocket) which makes it more difficult to extract the quasi- 
particle peak from the background. 

There are several issues which are not resolved in this 
paper and require further study. First, the photoemis- 
sion data indicate that the quasiparticle peak in the spec- 
tral function is broad already at optimal doping and be- 
comes even broader as the system moves towards half- 
filling. We have shown in Sec. II A that at exactly half- 



filling, the broadening of the quasiparticle dispersion is 
due to the interaction with short wavelength spin fluc- 
tuations (long wavelength fluctuations do not contribute 
because of the vanishing vertex). In our analysis of the 
spin-fermion model with a cutoff luq in the susceptibility, 
we completely neglected these fluctuations in the para- 
magnetic phase (we remind that in our discussion on 
the cancellation of the higher-order self-energy diagrams, 
we restricted with only the leading terms in the expan- 
sion in the bosonic frequency and the momentum shift 
from Q. We do not believe that the qualitative features 
of the Fermi-surface evolution will change if we include 
short-wavelength spin fluctuations, but they can be im- 
portant for a quantitative analysis. For example, the 



FLEX calculations, which do not yield precursors of the 
SDW state, also show that T c is reduced in underdoped 
cuprates simply because fermions become less coherent 
as one approaches half-filing f64|| . 

Another issue which requires further study has 
emerged from the photoemission data from the Stanford 
and Argonne groups |7]|^] which found that, in under- 
doped cuprates, A(u) for k near (0, n) not only has a 
broad peak at about 2 J, as it should be if the precursors 
of the SDW state are present, but also drops rapidly at 
frequencies of about SOmeV. There is no such drop for 
the data taken along the Brillouin zone diagonal. When 
the temperature is lowered below T c , the spectral func- 
tion acquires a narrow peak at exactly the same posi- 
tion where the drop has been observed |]65f . It is tempt- 
ing to associate this new feature with the precursors of 
the d— wave pairing state. Emery and Kivelson proposed 
a scenario for d— wave precursors in which proximity to 
the antiferromagnetic instability is not particularly rele- 
vant [[36| . A similar scenario has been suggested by Ran- 
deria et al. |57|. From our perspective, an important 
point is the observation by Shen et al. [ pl| that both, the 
high energy peak and the drop in A(uj) at low energies 
disappear at about the same doping concentration, i.e., 
the two features are likely to be correlated. This observa- 
tion poses the question whether precursors of the SDW 
state can give rise to precursors of the d— wave pairing 
state. At the moment, we do not know the answer to this 
interesting question. 

Finally, the reasons why Luttinger's theorem does not 
work above the critical value of the coupling are still not 
completely clear to us. Recently, we considered in de- 
tail why Luttinger's proof for / = J Gf u udH/du> = 
does not work in the magnetically-ordered state 68 1 . We 
found that though a formal application of Luttinger's ar- 
guments yields I = to all orders in the SDW gap A 
in perturbation theory, the frequency integrals contain 
a hidden linear divergence and have to be regularized. 
When this regularization is done, one obtains that / is 
finite, as it should be to recover Eq. (|Tl|) , and for small A 
behaves as I oc A 2 . We are currently studying whether 
the same reasoning can be applied above the topological 
transition in the paramagnetic phase. 

We conclude with a final remark. In this paper, we 
studied the Fermi surface evolution as a function of the 
coupling strength at zero temperature. There exist, how- 
ever, a number of experimental data which show temper- 
ature crossovers in various observables in the underdoped 
cuprates, most noticeably in the NMR relaxation rate, 
the uniform susceptibility and the resistivity. Recently, 
Pines, Stoikovich and one of us (A. Ch.) argued that 
these crossovers are related to the thermal evolution of 
the Fermi surface (or, more accurately, to the evolution 
of the quasiparticle peak in the spectral function) [f39| . 
We refer the interested reader to Ref. |39| for a detailed 
discussion of this issue. 
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APPENDIX A: THE SELF-ENERGY FOR 
ARBITRARY a k 

The explicit computation of Eq.(^8j) yields for arbi- 
trary a k = \v k+Q /c sw \ 



D 



+ - a l A o + %+Q I V 1 -^ 

\e k+Q \ + ^1-a 2 /e 2 +A 2 a 2 



(A5) 



if |efc+Q| < uj Q a k and 



B = 




V 1 - a W^o + ^o + + - a 2 ) 



(A6) 



2 \J (efc+Q ~ ium) 2 
~9 e ff — = : 

1 

(y/v$ + Al-A )y/l-ai 



if |efc+o| > LOQa k . We plotted this result (for a k — 1) in 
Fig. life. 



x A (Al) 



where 



and 



A = In 



A = c su ,/£, 9 ef f=g a^—^ y 

7T 



V(e fe +Q - iw m ) 2 + ^ 2 (1 - a 2 ) + V(l - «fe)(^o + A o) 



\/(efc+Q - «w m ) 2 + A ov /l - a; 



(A2) 



if £fc+Q + W^/o > a fc w 0> and 

^ = ln V (£fc+Q - »^m) 2 + £fc+Q + /q(l -~gfc) 

a (V (efe+Q ~ «^ m ) 2 + Aoy/l - d{\ 
| V^ 1 - a lWk+Q + ^ 2 » /q + « 2 cj 2 ) 
a k f V (efc+Q - iw m ) 2 + A ov /l - a£j 

if e fe+Q + w 2 J a > a\ujl with 

(l + a fe ) 2 -3 



(A3) 



fa 



a k 



For afc = 1 we recover the expressions presented in 
Eqs.© and ©. 

We also can extract th e qu asiparticle residue at the Fermi 
surface. Substituting (Al) into 



Z k 



= 1 



we we obtain 
1 



Z k 



1 



«9E 

duj 



lo 2 + A 2 - A ) (\e k+Q \ + A o yT 
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